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Abstract. Einstein's field equations for a spherically symmetric metric coupled 
to a massless scalar field are reduced to a system effectively of second order in 
time, in terms of the variables /i = m/r and y = (a/ra), where a, a, r and m 
are as in [W.M. Choptuik, "Universality and Scaling in Gravitational Collapse of 
Massless Scalar Field", Physical Review Letters 70 (1993), 9-12], Solutions for 
which 11 and y are time independent may arise either from scalar fields with ipt = 
or with (f>s = but rf> linear in t, called respectively the positive and negative 
branches having the Schwarzschild solution characterized by <^ = and /is + /i = 
in common. For the positive branch we obtain an exact solution which have 
been in fact obtained first in [I.Z. Fisher, "Scalar mesostatic field with regard for 
gravitational efll'ects" , Zh. Eksp. Tear. Fiz. 18 (1948), 636-640, gr-qc/9911008 
and rediscovered many times (see D. Grumiller, "Quantum dilaton gravity in 
two dimensions with matter", PhD thesis, Technische Universitdt, Wien (2001), 
gr-qc/0105078 I and we prove that the trivial solution = is a global attractor 
for the region fis + /i > 0, /i < 1/2. For the negative branch discussed first 
in [M. Wyman, "Static spherically symmetric scalar fields in general relativity". 
Physical Review D 24 (1981), 839-841] perturbatively, we prove that /i = is a 
saddle point for the linearized system, but the non- vacuum solution = 1/4 is a 
stable focus and a global attractor for the region /is + /J- > 0, fJ. < 1/2. 

Keywords: Choptuik spacetime; Spherically symmetric static solutions; Massless 
scalar fields 

PACS numbers: 83C15, 83C20 
1. Introduction 

The initial value problem for Einstein's equations with massless scalar field was studied 
analytically by Christodoulou in ^ and [21 where it was shown that "small initial 
data" disperses while for "large initial data" the end state is a black hole surrounded 
by vacuum. The work of Choptuik on the numerical search for "critical initial data" 
that would separate these two types of behavior led to the discovery of the "threshold 
phenomenon" . For various one parameter families of initial data, where p denotes the 
amplitude or the width of pulses of various shapes, it was observed that for p ^ 
the end state of the evolution is a flat space-time (external field dissipates), while for 
p ^ p, the end state of the evolution is a black hole. For p near p^,, where the black 
hole formation first starts, the black hole mass is proportional to {p — p*)"' where the 
constant 7 is nearly equal to 0.37 independently of the initial data. Furthermore, the 
critical solution exhibits a universal scaling behavior in the sense that it is periodical in 
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suitable time and space coordinates with a period nearly equal to 3.44 for all families. 
Threshold behavior was studied afterwards for a number of configurations and similar 
types of critical behavior have been observed. A detailed overview of the literature on 
the threshold phenomena can be found for example in 2]. 

The field equations used in numerical studies in the literature form a consistent 
but overdetermined system. In this paper we shall make use of the complete set of 
field equations as in jSj to decouple the equations for the metric functions from the 
equations for the scalar field. These equations form a normal system and they are 
effectively of second order in time. In logarithmic coordinates, they have no explicit 
dependency on the independent variables, as expected from the symmetry analysis 
of the original equations. The second order hyperbolic equation for the mass density 
involves a parameter that can be either plus or minus one, hence the decoupled system 
has positive and negative "branches" . Besides its suitability for perturbative studies, 
in the time independent case, this system reduces to an autonomous second order 
system and allows a rigorous analysis of the negative branch, which is not transparent 
in the original formulation. 

When the field equations are written in that form, a non-flat, special solution for 
the negative branch is obvious: fi = j, or m = jr, for which the external field is linear 
in time. Thus surprisingly, time varying external fields give rise to static metrics. A 
perturbative study of the full system is deferred to a future work. In Sections 4 and 
5 we restrict ourselves to the static case, which can be realized either as the ones 
corresponding to time invariant external fields or the ones with (ps = but 4>t 0. 
For the positive branch we obtain an exact analytic solution while for the negative 
branch, we prove that all solutions starting in the region /is + /x > 0, /x < 1/2, stay in 
that region and the constant solution /i = -j is a global attractor. Although there is 
no non- vacuum static black hole solution, for large values of (j)s oi (j)t, IJ- gets very close 
to 1/2 and decays rapidly to negative values, as depicted in Figure 3b. The trivial 
solution n = and the Schwarzschild solution with /i^ + /i = are of course common 
the both branches. 

It was pointed out by D. Grumiller in a private communication that the solution 
for the positive branch dates back to Fisher as discussed in detail in [7] and |S], 
Fisher's solution have been rediscovered in [HI-CSl ^^'^ errors of the original equations 
28-29 in |H] were corrected by Grumiller [Jj. In addition we may cite the works 
of Bronnikov and Shikin on the characterizations of the global properties of 
static, spherically symmetric configurations in various dimensions. As opposed to 
the popularity of the positive branch, the second class of static solutions that we call 
the negative branch is first noticed by Wyman where it is studied perturbatively. 

2. The Field Equations 

Let M be a four dimensional Lorentzian manifold and g be a spherically symmetric 
metric on M. In the coordinate system adopted in j^l, the local coordinates are 
xi^ = (t, r, 9, If) with 

< 9 < t:, < if <2'k, r > To > 0, t > to> 0. 

The metric is given by 



ds^ = 



-a^ ^2 ^^2 ^2 ^ ^2^-^20 ^^2 



(2.1) 
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where a = a{r,t), a ~ a{r,t). 

The Einstein's equations coupled to a scalar, static particle are 



where 



T^u = df,(t)d^(j) - -gp,u{d\(t)d 4>) 



and (j) satisfies the wave equation 

with g — dct((7^y) 16 . The complete set of field equations are given below, where we 
have written k = Stt and omitted an equation which was identically satisfied. 



a 



Or Oir C 

a a 

1 = Anr 

a a 



0, 



\a / 1 v a / r 



(2.2a) 

(2.26) 
(2.2c) 
(2.2d) 



Remark 2.1. The system above is over-determined but integrable in the sense that 
the equations for at and ar are compatible, as noted in The degree of freedom 
in the system appears to be three functions of a single variable, corresponding to the 
initial data for (jjt and a. This last freedom is actually redundant, since it can be 
eliminated by redefining t. The use of the complete system will allow the decoupling 
of the equations for a and a from the equation for (j). 

The literature on the threshold phenomenon works with the system of equations 
(2.2b-d), starting with initial data for (j). We shall use the function m{r,t) defined 
in terms of a(r, t) as our primary variable and introduce logarithmic coordinates to 
obtain a normal (i.e. not overdetermined) system with no explicit dependency on the 
independent variables. 

We start by defining the variables 



a 
a 



1 



2m 
r 



(2.3) 



Note that as r > 0, we need to impose r — 2m > to ensure the positivity of a^ . The 
equations (2.2a-d) can be rearranged in the form 



z 



2m 
r(r — 2m) ' 



r[r — 2m) 



47r0t(/)r, 



(2.4o) 
(2.46) 
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rrir 



27r 



r(r — 2m) 



(2.4c) 



(2.4rf) 



Finally passing to logarithmic coordinates we obtain an autonomous system, as given 
below in Proposition 2.1. 

Proposition 2.1. Let 

s = ln(r), z = ry, m = r^. (2.5) 
Then, the system of equations (2.2a-d) is equivalent to the system of equations 



^=-2- 

y 

1 Mt 
27r 1 - 2/Lt 



1 



1-2/x' 



27r 1 - 2m '^^ y2 ' 



— ) = 2y(/>s + {y<l>s)s, 



where Hs+ (J- o-nd Ht/y are constrained by 



Ms + M 
1 - 2n 



>0, 



1 {fis + m) ± ifit/y) 



27r 



1 - 2/1 



± — I > 0. 



(2.6a) 
(2.66) 
(2.6c) 
(2.6d) 

(2.7) 



Remark 2.2. In addition to the translation invariance, the system (2.6a-d) admits 

the symmetry 

Thus y can be scaled by an arbitrary function of t only, hence its time evolution is 
determined in terms of the time evolution of /x. 

Remark 2.3. A special case of (2.8) is the transformation 

{t,y)^ {T^\nt,y = ty). (2.9) 

Hence one can use (t, y) as well as {t,y) in (2.6a-d). If y is periodic in r, i.e., 
y(r,s) =y(T + T,s). Then 

e'^2/(e'^, s) = e^e'^2/(e'^e^, s), 

or 

y{t,s) = e^y{e^t,s), 
hence periodic solutions lead to the "echoes" of the original solution. 

The regions determined by the conditions = ±^ and + /l( = are shown 
in the Figures la and lb respectively. As the equations (2.6a-d) are invariant under 
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the discrete symmetry (pg — > —cpg it is sufficient to work with ^ in regions I and II for 
example. Also for 1 — 2/i > only region /* in Figure lb is allowed. As it will be seen 
in the next section, the half plane consisting of the regions / and // in the {4's>4't/y) 
plane is mapped to the quarter plane I* as a double cover. We shall call these as 
the "positive" and "negative" branches. The boundary of the regions / and // is 
characterized by 0s = ±0t/t/ or equivalently with + jj, = ±ij,t/y. The decoupled 
equation (3.1b) below admits non- vacuum solutions fj,s+IJ' = ^jJ-t / V, but when we take 
into account the wave equation for we get 0s = ^t/y = 0, hence /^s + M = = 0- 



3. Reduction of the field equations 

We shall eliminate from (2.6a-d) and obtain a system of equations for and y. The 
proof of the proposition below is essentially a straightforward computation. We have 
introduced the auxiliary variables X and Y below in order to make the proof more 
readable. 



Proposition 3. 1 . For 



(if )^ 7^ the system of equations (2.6a-d) is equivalent 



to the system of equations given below. 



y 



£1 = _2 -I 

y (i-2/x)' 



(3.1a) 



Ms 



y\y 



(1 - 2^) 



„2 

{lis + nf + ills 



±2e 



(/Zs + m)^-^ 

y 

1 A*. 



1/2 



i^l + {<t^t/yf = 
20s(0t/2/) = 



27r 1 - 2/x' 

}_JM/y_ 
2nl- 2ijl' 



(3.16) 
(3.1c) 
(3.1d) 



where e = sgn{l — 2/i). 
Proof. Let 



X = ^l + {cf>t/yf 

Y = 2M4>t/y) = 



1 Ats + M 

" 277 1- 2/i' 

1 Mi/y 



±^ is equivalent to VX^ - Y"^ = 0. 



27r 1 - 2/i 

It can be seen that either of the conditions 0s 

In the following we assume that X"^ — Y'^ ^0. Note that /it = hence y = can be 
achieved either by taking 0s = or 0t = 0. These time independent solutions will be 
discussed in Section 4. Assuming for example 0s ^ 0, we obtain 



(l>t/y 



Y 
20s 



Substituting this in X we have 



-x<i>l + - 



0. 
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From this we can solve 

ct>l = ^[X±Vx^^]. 

Writing 

A = X± _ y2 

we have 

As noted at the end of Section 2, the signatures of and (f>t are irrelevant, hence we 
may take the positive signature. 

After substituting these and multiplying with A^/^, the differential equation for 
(p reduces to 

1 



y 



1a 1a 

-A + -As 



l-2n 2 

Substituting A and taking derivatives one can see that the coefficients of and 
are proportional, in fact we have 

i[l ± X{X- - Yr"'][Y,ly -X.] = ^^{X± ^/^^) T 

It can be seen that 



Hence it follows that 



Y, - Xt/y = 

1-2/i 



Xs - Yt/y + ± 2^X2 - F2 = 

which is (3.1b). • 
For the rest of the paper we assume that /ig + /x > 0, hence we take e = 1 in the 
equation (3.1b). 

4. Static metrics: Exact solution for the positive branch 

In this section we shall study solutions for which the metric is independent of t, hence 
static. This will allow nevertheless a time dependency in the scalar field. Putting 
= 0, the system of equations (3.1a-d) reduces to 

^ = -2 + -i-, (4.1a) 

y 1 - 2/i 

2 

l^ss + Ms + X _ 2^ l^^^l^^ ^'"^ ^ ^^'^^ + m) = 0, (4.16) 
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We shall first give special solutions with fis = k/j,, then obtain the exact solution 
for the positive branch and finally study the solutions for the negative branch in 
Section 5. 

Special solutions with /Xg = fc/x. 

For the positive and negative branches wc obtain respectively, 
positive branch: {k + l){k + 2)fi = 0, 
negative branch: {k + 1) [8/^^ + {k — 2)/i] = 0. 
Exact solutions for those special cases are given below with equations (4.2-5), where 

2/o and /io arc integration constants. 

1. /z = is a special solution for both positive and negative branches for which 
4>s = 4>t =^ and y = yoe~^. The line element is the flat metric 

ds^ = -yl dt^ + dr'^ + d0'^ + r^sin^e dcp^. (4.2) 

2. /i= 1/4 is a solution for the negative branch corresponding to (j)^ = and 
(j)t/y = 1/2-^7?. For this case y is a constant. Note that such a solution intersects the 
m = r/2 line only at the origin. Thus its perturbations may lead to what is observed 
as black holes with infinitesimal mass. The metric (2.1) becomes 

ds^ = -2y>2 ^^2 ^ 2 dr^ + dO'^ + r'^sin^O dip^. (4.3) 
The only nonzero component of the Ricci tensor is Ru = 22/^ and the Ricci scalar is 

R = 

3. A: = — 1, or + /z = corresponding to constant m is a solution for positive and 
negative branches. Here = /XoC"'* and y = ~ 2/i|. This is the Schwarzschild 
solution, it is the only possibility for intersecting the /i = 1/2 line. For increasing s 
the solution curves tend to the origin. For negative they escape to infinity in the 
reverse direction. When /Zq is positive, it has the standard interpretation as the black 
hole mass. The metric (2.1) becomes 

-1 

dr'^ + de'^ + dip^. (4.4) 

4. fc = —2 or /Xs + 2^ = is a solution for the positive branch. Then /i = /ioC^^''' and 
y = yo I -j^ 1^^^ I 1 — 2/x 1^/^. For increasing s the curves tend to the origin. The part 
of the curve for positive n is not allowed. The part for negative /x escapes to infinity 
in the reverse direction. The metric (2.1) becomes 

-1 

dr^ + r^ d9^ + r^sin^O dip^. (4.5) 

The only nonzero component of the Ricci tensor is Rrr = ^a(2^^°-T-^) ^"^^ Ricci scalar 
isii=-%. 

Exact solution for the positive branch. 



2/Xo 


dt' + 


^ _ 2/Xo 


r 




r 



ds^ 



-yldt" 
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For the positive branch (jit = and we write = (/ig. The field equations reduce 

1 4/i- 1 



to 



^ = -2 + 
y 1-2IJ, 1-2/i' 

fXs = + 27r(l - 2/i)V'^ 
^ ^ l_ 

We will obtain an analytic solution for this system. 

Proposition 4-1- For = the solution of the system (4.6a-c) is 



(4.6o) 
(4.66) 

(4.6c) 



1 

^=2 



''^ 47r ^ 47r 



- = IV'I IV' -pI""^ + 

Do 



r 
r„ 



where tjj = <ps, r = e^ , yo and Tq are constants, 



16.)=-(1 + , 



where 



S = 



Vc^ + 167r 



(4.7) 

(4.8) 
(4.9) 

(4.10a) 

(4.106) 
(4.10c) 

(4.11) 



Proof. Let u = {1 — 2iJ.)tp. Taking the derivative of u with respect to ip and using 
(4.6b,c), we get 



This is a separable equation which can be solved as 



1 + ctp — AttiP"^ 



where c is an integration constant. From the definition of u we obtain /z as a function 
of tj) as in (4.7). Note that there are vertical asymptotes a,t ip =p and ip = —q, with p 
and q as in (4.10a,b). In the region — g < ^ < p, is less than 1/2 and has a unique 
maximum at ipM = c/Stt where 



1 c 

fiM = At(V'M) = /"(c/Stt) = - 



2 c2 + 167r 2 



(4.12) 
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For |c| > 4:^/^T, jjl crosses the level /U = 1/2 twice at ij)'^ given by 
A typical solution curve is given in Figure 2. Note that 



(4.13) 



^-=0 = 



1 



i.e., the parameter c is the slope of the curve /x(V') at V = or as s — > oo. As we may 

assume t/; > 0, we shall restrict in the range Q < < p. For c negative and zero, 
is negative, ranging from — oo to as V' changes from p to 0. For < c < 4-\/7r it has 
a maximum at Vm but lies below the /x = 1/4 line. As c gets large, the solution curve 
is nearly a broken line following closely the borders of the region bounded by /i = 1/2, 
tp = p and tp = —q. For the critical value of c = 4y^, /xm = 1/4 at ipM — ^/^^/tt and 
p = 2^{~^ + V^)- Typical solution curves are presented in Figure 3a. 
To solve for y we write 



dy 

y 



So 



i>-p 4> + q, 



where A = 1 and Bi and i?2 are as in (4.10c). By integrating we obtain y as in 
(4.8). It can be seen that dy/dip vanishes at tp^ for which i^{tp^) = 1/4. Hence for 
|c| < 4-^7?, y is a monotone increasing function of tp while for |c| > 4-^7?, y has two 
critical points. As it can be seen from Figure 4, for large c, the variation of y with tp 
is quite sharp. 

Finally, substituting the expression of fi in (4.6c), we obtain the equation 



ATTtplip"^ 



A Ci 



■)/) '^P — p V' + 9 



where A = — 1 and Ci and C2 are given by (4.10c). Integrating this and exponentiating 
we obtain (4.9). It can be seen that r = is a monotone decreasing function of V'- 
Parametric plots of /x, y and r in terms of ip, are shown in Figure 4. • 

From the remark above on the interpretation of c and a qualitative study of the 
solutions we see that the rate of change of the mass distribution with respect to ip as 
s — > 00 determines the solution completely. For large c, there is nearly a black hole 
formation at s* determined by ip{s*) = tpM = c/Stt, in the sense that, the solution 
curve nearly touches the fx = 1/2 line but moves away from it to go eventually to large 
negative values as s ^ — 00 or as r — > 0. 

We were unable to find an exact solution for the negative branch. For this case 
we shall study the equation (4.1b) as a dynamical system. To compare both branches 
and to analyze the efficiency of numerical algorithms, we also studied the positive 
branch of (4.1b) as a dynamical system in the plane. Writing u = fi^, we can 
express the positive branch of (4.1b) 



Ms 









-2>v — 2fi 



l-2/i 



(l/2 + 3l2U + 2/x2) 



(4.14) 
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The only critical point is (0, 0) and the linearized part at the origin is 







■ 1 ■ 






D 


s 


-2 -3 







(4.15) 



The origin is a node and /Xs + M = and /x^ + 2/x = are exact solutions. Typical 
solution curves shown in Figure 5 are sketched using the exact solution. Since tp 
is a monotone function of s, one can as well think of the phase plane curves as 
parameterized by tp and make parametric plots of yU and fig- It can be seen that 
curves in the second quadrant, lying between + ^ = and fXg + 2^ = correspond 
to c ranging from — oo to 0, Hs + = being associated with c = 0. The curves 
beyond /x^ + 2/i = cross the = line at a finite value of the parameter and reach 
the origin as s — > oo. The expressions of fj, and fXg obtained from the exact solution 
implies that the origin is a global attractor for all solutions starting in the region 
defined by /x^ + M > and fj, < 1/2. 



5. Static metrics: Phase plane analysis for the negative branch 

Writing u = ^s, we can express second order equation (4.1b) for ^ as a dynamical 
system 



Ms 









V 

2 (,,2 



1-2/^ 



(t/2 + i^iv + 2/i2) 



(5.1) 







■ 


1 " 






V 


s 


2 


1 




V 



where v = jig. There are two critical points, (0,0) and (1/4,0). The linearized part 
at the origin is 

(5.2) 

To study the local behavior at (0, 0) we note that the eigenvalues of the system above 
are —1 and 2 hence the origin is a saddle point. The corresponding eigenvectors are 
(1, —1) and (1, 2), hence the solutions curves approaching the origin along v = —fi in 
either direction are repelled and move away along u = 2/x. 
The linearized part at the point (1/4,0) is 



(5.3) 



The eigenvalues of the system above are — 1 ± iV3, hence the point (1/4, 0) is a stable 
focus. The local behavior near the points (0,0) and (1/4,0) will be crucial in the 
proof of the global behavior. We will now prove that (1/4, 0) is a global attractor for 
solutions in the open half plane /Xg + > 0, /x < 1/2. 

Proposition 5.1. Let be a solution of (5.1) and D be the region bounded by 

> Q and 11 < 1/2. // (^(0), z/(0)) belongs to D, then the solution curve {fi{s),v{s)) 
remains in D for all s and 

lim (/x(s),Ks)) = (1/4,0). 







■ 1 ■ 








s 


-4 -2 







Proof. Let V be the velocity vector of the solution curve. We can express V as 

2 



V = 



(i^ + 2/x)(z/ + 2/x- i 



(5.4) 
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As 1 — 2/i is positive on D, the direction of V is determined by the signatures of u, 
V + 2/i and + — 1/2. This leads to a subdivision of D into five regions as given 
below. 



Region 








IJ'S 




-Di: 




fi + u > 0, 2/i - 


f < 0. 


> 


< 




> 0, 


2/i + i/>0, 1- 




> 


> 


^3: 


> 0, 


1 - 4/x - 2zy > 0, 


H < 1/2. 


>0 


< 


D^: 


v<0, 


I - Afj, - 2iy > 0, 


/i < 1/2. 


< 


< 


D5: 


u<0, 


H + i'> 0, 1 - 


4/x-2i^<0, /i<5. 


< 


> 



In the region Di, the velocity vector of the solution curve lies in the fourth quadrant. 
By uniqueness it cannot intersect the + u = line. Furthermore as the origin is a 
saddle point the solution curve necessarily intersects the 1/ + 2/U = line and enters in 
the region D2- 

The velocity vector of the solution curves in the region D2 lie in the second 
quadrant, hence each solution curve passes to the region D3 for increasing s. 

In the region Ds, the velocity vector is in the fourth quadrant. Here we will show 
that the boundaries v + 2iJ. = 1/2 and /x = 1/2 are never hit. 

In order to show that solution curves do not hit the boundary v = — 2/i + 1/2, we 
prove that the slope of their velocity vector is larger than —2. For this we need the 
inequality 

{u + 2n){v + 2/t - 1/2) < v{l - 2fj,) 

to be satisfied. The equality in the relation above defines a hyperbola, the half line 
z/ + 2/i = 1/2. /i < 1/4 lies in the region where the inequality is satisfied. Hence 
away from the i/ = level, the velocity vector V is surely pointing away from the 
u + 2^= 1/2 line. On the other hand as (1/4, 0) is a focus, the solution curves cannot 
hit this point as gets close to zero, hence necessarily, cross the vertical line /i = 1/4. 

To prove that solution curves do not hit the line /U = 1/2 at a point (/U, with 
IJ. > 1/4, we need to show that the velocity vector V points below the line joining 
(/i, v) and (1/2,0), i.e., we need to prove that 

I '^s/i^s \>\ -m) I • 
Since (1 — 2/i) and jJ-s = v are positive in this region, we should have 

\vs\ (1 - 2/x) - 2v'^ > 0. 

Substituting for Us, we obtain 

4/ij/ + {2v + 2/i)(4/i - 1) > 0. 

For /I > 1/4 this inequality is surely satisfied and solution curves reach the /t axis 
without hitting the /t = 1/2 line. 

In the region D4, the velocity vector belongs to the third quadrant, hence the 
solution curves hit the z^ + 2/i = l/2 line and pass to the region D5. In the region D5 
the velocity vector belongs to the second quadrant. By uniqueness of solutions, they 
cannot cross the n + u = line and the local behavior guarantees that they stay away 
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from the origin and cross the fi axis. Again by uniqueness the solution curves spiral 
around the point (1/4,0) which is a stable focus. • 

The integral curves for the numerical solution of the dynamical system 5.1 in the 
domain /i^ + /i > and /i < 1/2 are given in Figure 7. Expect for the line /ig + /i = 0, 
all solution curves in this region spiral around the point (1/4,0) as s —^ oo. As one 
moves to right, the solution curves get very close to the line /.j, ~ 1/2 but decay quickly 
to the level /i — 1/4. The variations of y and (j)t/y with respect to s exhibit rapidly 
decaying oscillations as shown in Figure 8. 

The period of the linearized part of the dynamical system (5.1) at the point (-1,0) 
is ~ 3.62. For the nonlinear system we could observe only 3 periods of oscillations. 

Starting with the initial conditions (/z = 0, /i^ = /x°), ranging from 10~^ to 10^, we 
computed consecutive crossings of the integral curves across fi — l/A line at fi^ — cii, 
i — 1, 4 and we evaluated the periods T^, i = 1,2,3 as the distance s between two 
consecutive crossings. Representative values for a^'s and T^'s are listed below. 



10-^: 




0.0923 


0.0022 


0.0001 






3.5350 


3.6280 


3.6070 


= 1 : 


0-1,0,2,03 


0.4698 


0.0089 


0.0002 




Ti,T2,n 


3.3510 


3.6180 


3.6380 


: 10^ : 


ai,a2,03 


567.4564 


0.0837 


0.0020 




Ti,T2,n 


4.7040 


3.5430 


3.6280 



As /i^ goes to infinity the first crossing level oi goes to infinity i.e., lim^o^o^ ai = oo, 
while for i > 1 they tend to a limit that we denote a*, with 

lim Oi+i = a* = lim a^. 

For i = 1, 2, 3, we have a| ^ 0.092, Ss; 0.0022, = 0.0001. 

Although the first oscillation depends on the initial condition the second 
oscillation start near /is = 0.0022 — 0.0090 and have a quite stable period of 3.62, 
close to the period of the linearized system. 

Remark 5.1. Putting fitt — kj^s in the linearized part of the full system and writing 
/i = i + u we have 

Uss + (2 - k)us + 4m = 0. 

For A; = 2, the solution of linearized and nonlinear system are unrelated as predicted by 
Hartman-Grobman theorem For arbitrary < fc < 2 the linearized solution is of 
the form u — £"""6**"^ where a = 1 — | and w^+a^ = 4. As an interesting coincidence 
we observed that the period corresponding to A: = 0.3704 is T = j-ffgs — 3.44. 
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Figure 3a. Variation of /x with xjj = (pg for different values of c, for < V < P- 



Exact Static Solutions for Scalar Fields Coupled to Gravity in (3 + 1) -Dimensions 16 



0.5 







-0.5 



-1 



-1.5 



-2 















C.-4iF/ 




/ ^ r 






c=100 




12 



-2-101234 



Figure 3b. Variation of n with s for different values of c for — oo < s < oo. For 

c < 0, /i is monotone while for c > it has a unique maximum /iM at a finite s. For 
c large, yLt(s) makes a sharp peak at this critical point, as opposed to the nearly flat 
variation of /U with respect to <^s as seen in Figure 3a. 




Figure 4. The variations of \x, y and r with respect to ip = <ps for c = 6^/^T and 

< i/i < p. Recall that for c < 0, /i is monotone but has a unique maximum for 
c > 0. On the other hand, y is monotone for c < 4-^/7r but has two critical points for 
c > 4:^/^T while V' is monotone decreasing for all c. We have chosen the initial value of 
y as 2/o = 1- For increasing s, fi rises to its maximum /iM and decays to zero, while y 
starts from its initial value reaches its maximum, decays nearly to zero especially for 
large c and finally rises to oo very sharply. 
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Figure 5. The integral curves of the dynamical system (4.14) for the positive branch 

in the domain /ig + /i > and /i < 1/2. The lines fis + fJ- = and /is + 2)1 = arc 
exact solutions. The half line /Ug + /U = 0, /z < corresponds to the limit c ^ — oo. For 
the solution curves in between this half line and /Ug + 2/z = 0, c ranges from — oo to 
zero. For these solutions /i is negative, and monotone increasing and reach the origin 
as s ^ 00. For c positive, all solution curves reach their maximum value hm < 1/2 
and arrive at the origin as s — > oo. Note the sharp changes near the = 1/2 line for 
large c. 
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Figure 7. The integral curves of the dynamical system (5.1) for the negative branch 

in the domain /is + /i > and fi < 1/2. The line /i^ + /i = is an exact solution. 
All solution curves in this region spiral around the point (1/4,0) as s ^ oo. The 
curves sketched here are determined by their initial values (/i(0), /Lts(O)), which are 
respectively (-1,1.001), (-0.6,0.9), (0.2,1), (0.4,1), (0.49,1) respectively from right 
to left. Note that as one moves to right, the solution curves get very close to the line 
/i = 1/2 but decay quickly to the level /U = 1/4 making small amplitude oscillations 
around this level. 
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Figure 8. The variations of y, jj, and (pt/y with respect to s for the dynamical system 

(5.1) and Eqn.(2.6a). The curves sketched here are determined by their initial value 
{y{0) , IJ,{0) , IJ.s{0)) , which is (1,-0.6,0.9). Note that the solution of 4>t/y are obtained 
from Eqn.(2.6c) for cp^ = 0. 



